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Abstract 

A  simple  theoretical  model  is  presented  to  simulate  the  galvanostatic  discharge  behavior  of  the  Ni-composite  graphite  electrode.  The 
discharge  profiles  predicted  by  using  a  constant  diffusion  coefficient  (CDC)  and  by  a  varied  diffusion  coefficient  (VDC)  are  compared  in 
this  paper.  The  results  show  that,  the  VDC  model  can  be  simplified  to  the  CDC  for  discharge  rates  less  than  2C  for  a  5  pm  particle.  Also,  an 
approximate  analytical  solution  is  presented  for  VDC  model,  which  is  found  to  be  valid  for  discharge  rates  up  to  6C.  Exchange  current  and 
diffusion  coefficient  for  the  lithium-diffusion  are  predicted.  ©  2001  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

Li-ion  batteries  are  now  the  preferred  power  sources  for 
cellular  phones,  portable  computers,  and  camcorders.  The 
distinguished  feature  of  lithium-ion  batteries  is  the  use  of 
intercalation  materials  for  both  positive  and  negative  elec¬ 
trodes.  Because  the  intercalated  species  must  transport  from 
deep  within  the  electrode  to  the  surface  to  participate  in 
electrochemical  reactions  that  ultimately  deliver  energies, 
the  particle  size  of  the  intercalation  electrode  plays  an 
important  role  in  determining  the  energy  performance  of 
the  Li-ion  cell  [1].  For  example,  particles  too  large  give  rise 
to  low  energy  densities  at  even  moderate  discharge  rates 
[1,2].  In  contrast,  particles  too  small  tend  to  be  difficult  to 
incorporate  into  electrode  manufacturing  processes.  In  the 
case  of  graphite  negative  electrode,  particles  too  small  also 
yield  a  high  surface  area,  which  eventually  leads  to  a  large 
irreversible  capacity  loss  arising  from  the  electrolyte  decom¬ 
position  [3,4].  Therefore,  the  particle  size  of  the  host  mate¬ 
rial  needs  to  be  carefully  designed  to  obtain  the  maximum 
battery  energy.  In  designing  lithium  rechargeable  batteries, 
some  other  parameters  such  as  the  electrode  porosity,  the 
lithium  diffusion  coefficient  in  the  solid  electrode,  and  the 
exchange  current  density  of  the  electrode  may  also  signifi¬ 
cantly  influence  the  battery  performance.  Hence,  a  mathe¬ 
matical  model  taking  into  account  these  material  properties 
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is  essentially  required  to  give  the  insights  of  the  physical 
processes  occurring  within  the  intercalation  host  and  to 
predict  theoretically  the  behavior  of  the  intercalation  host. 
This  in  turn  can  lead  to  a  better  development  and  design  of 
the  intercalation  materials. 

A  number  of  models  have  been  derived  for  analyzing  the 
behaviors  of  insertion  electrodes  [5].  West  et  al.  [6]  derived  a 
mathematical  model  that  accounts  for  the  coupled  transport 
in  the  electrode  and  electrolyte  phases  during  the  charge  and 
discharge  of  the  TiS2  porous  electrode.  Mao  and  White  [7] 
extended  West’s  model  by  including  the  mass  transport 
across  a  separator.  A  complete  model,  in  which  the  con¬ 
centrated  solution  theory  was  incorporated,  was  derived  by 
Doyle  et  al.  [8]  for  simulating  the  discharge  of  a  lithium 
electrode/polymer  electrolyte/TiS2  insertion  cathode  cell. 
Fuller  et  al.  [9,10]  and  Doyle  et  al.  [11]  presented  similar 
mathematical  models  for  simulating  the  galvanostatic  dis¬ 
charge  of  petroleum  coke/electrolyte/LivMn204  dual  lithium 
ion  insertion  cell,  for  studying  the  relaxation  phenomena  in 
petroleum  eoke/electrolyte/LivCo04  system,  and  for  exam¬ 
ining  the  behavior  of  the  plastic  lithium  ion  cells.  A  key 
approximation  used  in  all  these  models  is  that  the  diffusion 
coefficient  of  insertion  species  is  assumed  to  be  constant 
over  the  entire  concentration  range,  although  the  actual 
value  of  diffusion  coefficient  has  been  reported  to  be  con¬ 
centration-dependent  [12-15].  Recently,  Verbrugge  and 
Koch  [16,17]  incorporated  the  variation  of  diffusion  coeffi¬ 
cient  with  concentration  into  a  mathematical  model  to 
simulate  the  cyclic  voltammograms  and  the  galvanostatic 
intermittent  titration  process  of  the  lithium  intercalation  in 


0378-7753/01/$  -  see  front  matter  ©  2001  Elsevier  Science  B.V.  All  rights  reserved. 
PII:  S0378-7753(00)00657-l 


V.R.  Submmanian  et  al.  /  Journal  of  Power  Sources  96  (2001)  396—405 


397 


Nomenclature 

As 

surface  area  of  the  electrode  per  unit  volume 
(cm2/cm3) 

initial  (reference)  lithium  concentration  in 
carbon  insertion  host  (mol/cm3) 

De 

D\ 

diameter  of  the  electrode  (cm) 

chemical  diffusion  coefficient  of  lithium  in 

solid  graphite  at  the  infinite  dilute  state  (cm2/ 
s).  Found  out  to  be  (>  1.25  x  10“9  cm2/s) 

F 

Faraday’s  constant  (96487  C/equivalent) 

f(X i) 

function  of  diffusion  coefficient  dependence 
on  the  concentration 

h 

local  current  density  on  the  particle  surface 
(A/cm2) 

^  o.  re  1’ 

local  reference  exchange  current  density 
(A/cm2) 

/ 

macro  current  density  of  the  electrode  (A/g) 

Axref 

macro  reference  exchange  current  density 
(A/g) 

j 

integer  number  in  Eq.  (13) 

k 

integer  number  in  Eq.  (4) 

L 

thickness  of  the  electrode  (cm) 

me 

mass  of  the  electrode  (g) 

Mw 

molecular  weight  of  graphite  (12  g/mol) 

N 

integer  number  in  Eq.  (4),  (here  N  =  8) 

Q 

amount  of  charge  passed  to  the  graphite 
electrode,  Q  =  It  (Ah/g) 

r 

radial  coordinate  of  a  spherical  particle  (m) 

R 

universal  gas  constant  (8.3143  J/mol/K) 

Rp 

particle  radius  (cm) 

t 

time  (s) 

u 

electrode  potential  relative  to  Li+/Li  (V) 

ve 

geometric  volume  of  the  electrode  (cm3) 

X 

stoichiometric  number  in  LiTC6 

X, 

fractional  occupancy  of  lithium  in  graphite 
host  (0  <  X\  <  1) 

n 

initial  fractional  occupancy  of  lithium  ion  in 
graphite,  obtained  by  the  ratio  of  the  full 
discharge  capacity  to  the  theoretical  capacity 
of  graphite  (372  mA  h/g) 

Ms 

fractional  occupancy  of  lithium  ion  on  the 
particle  surface 

J 

dimensionless  radial  coordinate 

Greek  letters 

P 

symmetric  factor  (=  0.5) 

£ 

porosity  of  the  electrode 

7i 

activity  coefficient 

h 

root  of  tan  a  =  a.  j  =  1 , 2, . . . ,  oo 

P 

specific  density  of  graphite  (g/cm3) 

T 

dimensionless  time 

rk 

polynomial  fitting  parameter,  k  =  1,  2, . . . ,  8 

Qk 

self-interaction  coefficient,  k  =  2, 3, . . . ,  N  (in 
this  work  N  =  8) 

<P  open-circuit  potential  relative  to  Li+/Li  (V) 

(p°  concentration-independent  potential  of  the 

graphite  electrode  relative  to  Li+/Li  (V) 

0S  thermodynamic  potential,  defined  in  Eq.  (4) 
(V) 

¥  dimensionless  parameter  denoting  the  ratio  of 

diffusional  resistance  of  lithium  ions  in  the 
particle  to  the  interfacial  kinetic  resistance, 
defined  in  Eq.  (15) 


the  carbon  fiber,  however,  how  the  varied  diffusion  coeffi¬ 
cient  affects  the  predicted  result  was  not  illustrated  in  their 
model. 

In  this  work,  a  simple  mathematical  model  for  the  galva- 
nostatic  discharge  process  of  lithium  in  graphite  porous 
electrode  is  developed.  The  model  developed  is  further  used 
to  analyze  the  effects  of  design  parameters  on  the  lithium 
intercalation  into  graphite. 

2.  Experimental 

The  model  simulation  was  performed  on  the  system 
of  a  porous  graphite  pellet  electrode,  which  was  used  as 
a  lithium  intercalation  host.  The  porous  electrode  was 
made  of  Ni-composite  graphite  powder,  which  was  synthe¬ 
sized  by  electroless  deposition  of  10  wt.%  thin  nickel-com¬ 
posite  coating  (approximate  1  nm  in  thickness)  [18]  on 
the  surface  of  graphite  KS10  (10  pm,  Timcal  America). 
The  Ni-composite  coating  on  graphite  was  utilized  to  re¬ 
duce  the  irreversible  capacity  loss  of  graphite  in  the  initial 
cycle.  Since  the  coating  does  not  impede  the  lithium  inter¬ 
calation  into  graphite,  the  presence  of  Ni-composite  coating 
does  not  affect  the  predicted  charge  performance  of  the 
graphite. 

The  graphite  pellet  electrode  was  prepared  by  mixing  the 
Ni-composite  graphite  powder  with  6  wt.%  poly(vinylidene 
fluoride)  (PVDF,  Aldrich)  powder  and  l-methyl-2-pyrroli- 
dinone  solvent  on  a  stainless  steel  current  collector  and  dried 
under  vacuum  at  150°C  for  12  h.  The  obtained  pellet  is  about 
9.9  mg  in  weight,  1.26  cm  in  diameter  and  94  pm  in  thick¬ 
ness.  The  volume  of  this  pellet  electrode  is  calculated  to  be 
0.0117  cm3.  The  electrolyte  consists  of  1  M  LiPF6  in  a  1:1:3 
mixture  of  PD/EC/DMC  with  less  than  15  ppm  H20  and 
80  ppm  HF  (EM  Inc.). 

The  galvanostatic  charge-discharge  experiment  was  con¬ 
ducted  in  a  Swagelok  three-electrode  cell,  in  which  the 
graphite  pellet  served  as  a  working  electrode  and  lithium 
metal  foil  was  used  as  a  reference  electrode  as  well  as  a 
counter  electrode.  The  experiments  were  carried  out  at  25°C 
using  the  EG&G  PAR  Potentiostat/galvanostat  Model  273 A 
driven  by  the  CorrWare  software  system  from  Scribner 
Associates,  Inc.  The  galvanostatic  charge-discharge  data 
were  collected  after  cycling  three  times  to  ensure  that  stable 
and  reproducible  behavior  had  been  obtained. 
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3.  Diffusion  model 


The  Ni-composite  graphite  particles  in  the  porous  gra¬ 
phite  pellet  electrode  are  assumed  to  be  in  spherical  shape. 
The  electrochemical  reaction  occurring  at  the  individual 
particle  surface  is  described  as 

Charge 

xLi++xe“+C6  ^  LixC6  (1) 

Discharge 

where  x  is  a  stoichiometric  number  of  lithium  intercalated 
into  the  graphite  host.  Upon  charging  the  graphite,  lithium 
ions  diffuse  from  the  electrode  surface  into  the  graphite 
insertion  site.  The  reverse  process  occurs  during  the  dis¬ 
charging  process.  To  simplify  the  analysis,  only  the  electro¬ 
chemical  intercalation/deintercalation  reactions  shown  in 
reaction  1  are  considered  in  this  model.  The  side  reactions 
such  as  the  electrolyte  decomposition  and  lithium  metal 
deposition  are  neglected.  The  diffusion  of  lithium  ions  inside 
the  spherical  particle,  can  be  represented  by  Fick’s  second 
law: 


dXi 

~di 


D\  d 
r 2  dr 


(2) 


where  X{  denotes  the  fractional  occupancy  of  lithium  ions  in 
total  available  sites  of  the  insertion  host,  X\  =  <2/372,  where 
Q  and  372  mA  h/g  are  the  actual  capacity  and  the  theoretical 
capacity  of  lithium  in  graphite,  respectively.  The  value  of  A, 
ranges  from  0  to  1.  X\  =  1  refers  to  the  state  that  lithium 
occupies  all  the  available  sites,  which  corresponds  to  the 
capacity  of  372  mA  h/g  and  r  is  the  radial  distance  from  the 
center  of  the  particle  and  f{X{)  is  the  concentration  depen¬ 
dence  of  the  diffusion  coefficient. 


The  open  circuit  potential  (<P)  of  the  lithiated  graphite 
electrode  relative  to  the  metallic  lithium  can  be  formulated 
as  a  function  of  the  lithium  fractional  occupancy  at  the 
particle  surface  [16,17]: 


<P 


1  —  Xi.s 


(3) 


where  &°  refers  to  the  electrode  potential  at  the  infinite  di¬ 
lute  state  of  lithium  in  graphite  [16],  i.e.  at  Aj  s  — >  0  and  R,  T 
are  the  universal  gas  constant  and  temperature,  respectively. 

The  thermodynamics  of  lithium  intercalations  into  a  solid 
insertion  host  departs  from  the  Nernst  equilibrium  because 
the  concentrations  of  the  intercalated  species  are  higher  than 
the  dilute  value  assumed  at  the  standard  state  [10,16].  The 


parameter  Qk  in  Eq.  (4)  is  determined  by  fitting  Eq.  (3)  to  the 
experimental  curve  of  the  open-circuit  potential  of  graphite 
versus  the  lithium  occupancy  using  a  least  square  polyno¬ 
mial  regression  method.  Rewriting  Eq.  (3) 


<P  =  <PS+  —  In- — — ,  where  <PS  =  <P°  -  S^—kX^ 
F  Xls  f-jF  "■ 


(4) 

Since,  there  is  a  difficulty  in  collecting  the  actual  equili¬ 
brium  data  of  the  carbon  electrode  in  a  finite  period  of  time, 
the  charging  potential  profile  at  a  low  charge  current  has 
been  assumed  to  be  an  equilibrium  potential  profile  and  used 
to  determine  the  thermodynamic  parameters  for  the  model¬ 
ing  simulation  [8-11].  The  open-circuit  potential  profile  is 
obtained  by  charging  the  graphite  system  at  C/40  rate 
(9.3  mA/g),  which  is  lower  than  the  rate  of  C/24  (i.e. 
15.5  mA/g)  known  to  be  sufficiently  slow  to  attain  the 
near-equilibrium  condition  on  graphite  [1].  Fig.  1  presents 


0.0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9 

x  in  LixC6 


Fig.  1.  Open-circuit  potential  as  a  function  of  the  composition  (jc)  in  Li  ,C6. 
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Table  1 


Activity  coefficients 


Material 

Our  result 

Graphite  powder 

Verbrugge’s  result 

Partially  graphitic  carbon  fiber 

0°  (V) 

1.3935 

0.817 

02/F  (V) 

10.6885 

0.9926 

03/F  (V) 

50.8584 

0.8981 

0fF  (V) 

148.8472 

-5.630 

05/F  (V) 

270.356 

8.595 

0JF  (V) 

296.0333 

-5.784 

07/F  (V) 

-178.537 

1.468 

0S/F  (V) 

45.4660 

the  experimentally  thermodynamic  data  of  graphite.  The 
thermodynamic  data  for  partially  graphitic  carbon  fiber  by 
Verbrugge  [16]  is  also  plotted  for  comparison.  Note  that,  a 
plateau  at  the  potential  of  0. 1  V  versus  Li+/Li  which  char¬ 
acterizes  the  stage  phenomena  of  lithium  graphite  intercala¬ 
tion  compound  (Li-GIC),  is  observed  in  the  solid  dot  curve 
for  graphite,  while  no  significant  plateau  displays  for  par¬ 
tially  graphitic  carbon  fiber.  The  solid  dot  data  of  graphite 
are  fitted  to  Eq.  (4)  using  a  polynomial  regression  method 
and  the  resulted  fitting  curve  is  depicted  as  a  solid  line  in 
Fig.  1.  The  values  of  potential  <P°  and  the  interaction 
coefficients  Qk/F  (k  =  2-8)  are  determined  and  summarized 
in  Table  1.  The  values  obtained  by  Verbrugge  [16]  are  also 
included  in  Table  1  for  comparison.  The  differences  in  the 
potential  and  interaction  coefficients  probably  arise  from  the 
different  materials  used  for  lithium  intercalation.  As  seen  in 
Fig.  1,  the  lithium  intercalation  in  graphite  powder  electrode 
mainly  occurs  at  0.1  V  versus  Fi+/Fi  while  the  lithium 


intercalation  in  partially  graphitic  carbon  occurs  from 
0.8  V  versus  Fi+/Fi  and  continues  over  the  potential  range 
of  0. 1-0.8  V  versus  Fi+/Fi.  These  different  intercalation 
potentials  could  result  in  the  different  values  of  (p°  and  QkIF 
for  the  partially  graphitic  carbon  and  the  graphite  powder  as 
shown  in  Table  1. 

The  function  of  diffusion  coefficient  dependence  on  the 
concentration  is  obtained  [17]  as 


m i) 


_  DijXi)  _ 

D\ 


dlnqq 
din  Vi 


-Xi(I-Xi) 


F  d& 
RTdXi 


(5) 


where  yT  is  the  activity  coefficient  of  lithium  in  the  solid 
graphite  host.  The  thermodynamic  data  of  Ni-composite 
graphite  in  Fig.  1  were  differentiated  with  X\  using  a  five- 
point  finite  difference  method  (accurate  to  the  order  h4)  to 
obtain  the  values  of  d<P/dXi  at  different  concentration.  The 
values  of  /fV,)  at  various  Xj  were  calculated  from  Eq.  (5)  and 
plotted  as  a  solid  line  in  Fig.  2.  Note  that,  three  peaks  in  Fig.  3 
are  observed  at  the  concentration  Vj  of  about  0.09,  0.22,  and 
0.51.  These  Vj  values  as  shown  in  Fig.  1  are  in  the  ranges  that 
the  potential  quickly  drops,  which  corresponds  to  the  stage 
transform  of  lithium  intercalates.  This  result  indicates  that 
the  significant  variation  of  the  lithium  diffusion  coefficient 
in  graphite  occurs  at  the  state  of  the  stage  transform.  The 
solid  line  of  Fig.  2  is  separated  into  five  pieces  and  is 
regressed  using  a  polynomial  regression  method  to  obtain 
the  expression  of  /fV,) 


/(X,)  =  ( 


k=l 


(6) 


Fig.  2.  Diffusion  coefficient  as  a  function  of  fractional  occupancy. 
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Q  (mAh/g) 

Fig.  3.  Comparison  of  the  experimental  and  the  simulated  discharge  profile  at  the  discharge  rate  of  C/8  (i.e.  46.5  mA/g). 


The  fitted  curve  is  depicted  as  the  dotted  line  in  Fig.  3.  The 
corresponding  fitting  parameters  of  r k  are  summarized  in 
Table  2. 

The  discharge  process  starts  from  the  fully  charged  state 
in  which  all  the  lithium  ions  uniformly  distribute  in  the 
graphite  particle  with  a  fractional  occupancy  of  X{J.  The 
value  of  Xj*  is  obtained  by  the  ratio  of  the  full  charge  capacity 
to  the  theoretical  capacity. 

at  t  =  0,  Xi  =  Xj°  (7) 


and  the  boundary  conditions  are 
dX] 

at  r  =  0,  =  0 

or 

n  n  C)X I  lr\ 

at  r  =  D?C°  — ! -  = - 

p’  1  s  dr  Ff(Xi) 


(8) 

(9) 


where  F  is  the  Faraday’s  constant.  C°  is  the  initial  (refer¬ 
ence)  lithium  concentration  in  graphite  host,  calculated  with 


the  fractional  occupancy  of  lithium  at  fully  insertion  state 
(X[J),  the  specific  density  (p)  and  the  molecular  weight  (Mw) 
of  graphite  from  the  expression  [19] 


C°  = 


Xjp 

6  Mw 


(10) 


where,  6  refers  to  the  stoichimetric  number  of  carbon  in 
Li^Cg. 

In  Eq.  (9),  ip  denotes  the  local  current  density  on  the 
particle  surface  in  terms  of  A/cm2.  The  value  of  ip  is 
calculated  with  respect  to  the  macro  current  density  of 
the  electrode,  I  (A/g),  the  mass  of  the  electrode,  me  (g), 
the  surface  area  of  the  electrode  per  unit  volume,  As  (cm2/ 
cm3),  the  geometric  volume  of  the  electrode,  Ve  (cm3),  and 
the  electrode  porosity  e,  [20]  as 


.  _  /me 

'p  _  AsFe(l  -e) 


(11) 


Table  2 

Fitted  values  of  rk 


0<XI<  0.303 

0.303  <  X1  <  0.430 

0.431  <X1  <  0.530 

0.530  <  Xj  <  0.755 

0.775  <  XI  <  0.877 

r0 

0.773639563 

0.19 

247332.207972 

0.13 

-17977.4156 

r  i 

408.1776515 

0 

-1660960.478778 

0 

78323.91484 

r2 

-20856.54417 

0 

212827.412501 

0 

-109505.6582 

r3 

544047.5498 

0 

21964243.996622 

0 

51468.93968 

r4 

-7112522.79 

0 

-43752723.673841 

0 

-55839.23733 

r5 

49518548.04 

0 

-75525635.243337 

0 

130202.4323 

r6 

-187970332.8 

0 

367915176.568545 

0 

-58880.87905 

r7 

368041823.4 

0 

-472426675.845534 

0 

-55749.27299 

rs 

-291326139.5 

0 

210250890.426413 

0 

38051.47385 
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The  surface  area  of  the  electrode  per  unit  volume  (As)  is 
related  to  the  radius  of  the  spherical  particle,  Rp,  and  the 
electrode  porosity 


512  MB  RAM,  dual  Pentium  processor  is  less  than  1  min. 
The  same  program  can  be  used  to  obtain  the  discharge 
profiles  for  CDC  models  also. 


As 


3(1  —  e) 

Kp 


Three  types  of  models  are  described  here. 


3.1.  Constant  diffusion  coefficient  (CDC) 


(12) 


The  diffusion  Eq.  (2),  with  the  initial  and  boundary 
conditions  given  by  Eqs.  (7) — (9),  for  constant  diffusion 
coefficient  (CDC)  can  be  conveniently  solved  by  means 
of  variables  separation  [20-22]  as 


ipRp 

fd°ycs 


^py sin (Air/R p)  -£(tD0JRi) 

r  h  sin 


(13) 


where  the  f  is  the  root  of  tan  /.  =  /,  (/'  =  1.2....,  oo) 
[20-22].  Eq.  (13)  can  be  written  in  the  dimensionless  form 


=  X°  -  ¥ 


3t  +  -  (5y2  -  3)  -  -  V  f(4v)  e-^ 
10  ;  yUtfsin  (A,) 


(14) 


with  the  dimensionless  variables 


¥ 


ipRp 

FDfCs' 


tD\ 


(15) 


¥  denotes  the  ratio  of  the  diffusional  resistance  of  lithium 
ions  in  the  particle  to  the  interfacial  kinetic  resistance.  When 
¥  <§  1,  the  diffusion  resistance  within  the  graphite  particle  is 
less  important.  On  the  other  hand,  when  ¥  >  1,  it  is  possible 
for  the  charge-discharge  of  the  electrode  to  be  limited  by  the 
diffusion  of  lithium  within  the  graphite  insertion  electrode. 
This  limitation  often  occurs  at  the  high  rate  of  early  cells 
using  nonporous  thin  film  [23].  The  lithium  fractional 
occupancy  at  the  particle  surface,  Xi  s,  is  obtained  by  sub¬ 
stituting  y  =  1  into  Eq.  (14) 


Xi,s  =  X°-¥ 


3t  +  0.2  — 


OO  I 

7=1  Aj 


-XU 


(16) 


3.2.  Variable  diffusion  coefficient  (VDC) 

When  the  diffusion  coefficient  varies  at  different  concen¬ 
tration  as  D\  =  Dj/(X i),  the  diffusion  Eq.  (2)  with  the  initial/ 
boundary  conditions  (Eqs.  (7) — (9))  is  solved  numerically 
using  finite  differences  in  y  and  integrating  the  resulting 
nonlinear  Ordinary  differential  equations  (ODEs)  using 
Maple  [24].  The  Maple  program  uses  40  node  points  in 
y-axis.  The  time  taken  for  this  program  in  an  833  MHz, 


3.3.  Variable  diffusion  coefficient-approximate  solution 
(VDCA) 


An  approximate  solution  for  this  concentration  dependent 
diffusion  case,  [25,26]  can  be  obtained  by  assuming  that  the 
concentration  inside  the  particle  is  parabolic  in  nature.  The 
details  regarding  the  derivation  can  be  found  in  [25,26]. 

.  1  5/(Xi,s)(l-A',,s)-y 

t  =  l5 - fxf) -  (17) 

The  accuracy  of  this  expression  depends  on  both  ¥  and 
f(X x).  Note  that  AI  S  is  given  by  a  non-linear  implicit 
algebraic  expression  in  Eq.  (17).  For  higher  values  of  ¥, 
one  can  use  higher  order  approximations  [25,26]. 


3.3.1.  Electrokinetics 

The  kinetic  of  charge-transfer  reaction  occurring  on  the 
particle  surface  is  described  by  the  Butler- Volmer  equation, 
which  relates  the  local  current  density  (L)  to  the  surface 
overpotential  (pf)  [17]  as 


/ 

ru -/>)*■„  i 

<  exp 

RT 

exp 


PF 

- 1 

RT 


(18) 


where  jl  is  the  symmetric  factor  ((>  =0.5  here)  and  /G  ref 
is  the  local  reference  exchange  current  density  in  terms  of 
A/cm2.  The  overpotential  (t]s)  is  specified  from  Eq.  (4) 


>h  =  U -<P 


U  — 


1-Xy; 
Xi, 


(19) 


Eqs.  (16)  and  (19)  are  the  governing  equations  for  the 
dependent  variables  XI  S  and  U  at  the  different  time,  when 
a  constant  diffusion  coefficient  (CDC)  is  assumed.  Upon 
discharging  at  a  time  t,  the  lithium  fractional  occupancy  at 
the  surface  of  particle,  XI  S,  is  determined  first.  Then,  from 
Eq.  (19)  the  electrode  potential,  U  is  obtained.  The  potential 
U  is  plotted  as  a  function  of  the  discharge  capacity  given  by 
Q  =  It. 

In  the  case  of  a  varied  diffusion  coefficient  (VDC),  the 
surface  concentration  distribution  with  time  is  obtained 
numerically  and  the  potential,  U  is  obtained  later  from 
Eq.  (19).  Alternatively  Eq.  (17)  is  used  for  finding  the 
surface  concentration  in  the  VDCA  model. 


4.  Results  and  discussion 

The  basic  parameter  data  used  for  the  modeling  simula¬ 
tion  are  listed  in  Table  3.  The  parameters  measured 
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experimentally  are  based  on  Ni-composite  coated  KS10 
graphite  and  the  manufactured  data  are  obtained  only  for 
bare  KS10  graphite.  The  simulation  process  stalls  from  the 
fully  charged  state  in  which  the  capacity  of  lithium  ions  in 
the  graphite  host  is  obtained  to  be  326  mA  h/g  at  a  charging 
rate  of  C/40,  i.e.  9.3  mA/g  (C  rate  refers  to  372  mA/g  of  a 
current  density).  The  initial  (reference)  lithium  fractional 
occupancy  (Aj’)  is  determined  to  be  0.877  by  taking  the  ratio 
of  this  capacity  over  the  theoretical  capacity  of  graphite 
(372  mA  h/g).  The  lithium  reference  concentration,  C°  is 


determined  to  be  0.027  by  applying  the  values  of  0.877  for 
A{3,  2.21  g/cm3  for  p,  and  12  g/mol  for  Mw  into  Eq.  (10).  The 
initial  guesses  for  the  diffusion  coefficient,  Dj,  and  the 
reference  exchange  current  density,  /o  re  r,  are  arbitrary. 
These  two  values  are,  then  predicted  accurately  using  the 
model  simulation. 

The  predicted  discharge  potential  profile  and  the  experi¬ 
mental  curve  as  a  function  of  discharge  capacity  are  plotted 
in  Fig.  3  for  C/8  rate.  The  dotted  line  is  experimental  data 
and  the  solid  line  is  the  simulated  result.  All  of  CDC,  VDC 


(b)  Capacity  (mAh/g) 


Fig.  4.  Comparison  of  the  predicted  discharge  profiles  by  the  CDC  model  and  the  VDC  model  at  different  discharge  rates. 
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Table  3 


Summary  of  the  parameters  used  in  the  simulation 


Parameter 

Value 

Resource 

Particle  radius  (pm) 

5 

Timcal  America  (manufacturer) 

Diameter  of  the  graphite  electrode,  De  (cm) 

1.26 

Measured 

Thickness  of  the  graphite  electrode,  L  (cm) 

9.4  x  1(T3 

Measured 

Mass  of  the  graphite  electrode,  me  (g) 

0.0099 

Measured 

Volume  of  the  graphite  electrode,  Ve  (cm3) 

1.17  x  10-2 

Measured 

Porosity  of  the  graphite  electrode,  (s) 

0.35 

Assumed 

Surface  area  per  unit  volume,  As  (cm2/cm3) 

3.9  x  103 

From  Eq.  (12) 

Specific  density  of  graphite,  p  (g/cm3) 

2.21 

Timcal  America  (manufacturer) 

Molecular  weight  of  carbon,  Mw  (g/mol) 

12 

Assumed 

Initial  lithium  fractional  occupancy,  Vj1 

0.877 

Measured 

Site  concentration,  Cs  (mol/cm3) 

0.027 

From  Eq.  (7) 

Symmetric  factor,  /I 

0.5 

[16] 

Diffusion  coefficient,  D]1  (cm2/s) 

io-“ 

Initial  guess 

Reference  exchange  current  density,  7ojef  (mA/g) 

10 

Initial  guess 

and  VDCA  models  give  the  same  profile  and  predict  the 
same  value  of  diffusion  coefficient  and  exchange  current. 
The  experimental  data  were  taken  every  12  min  from  a  cell 
that  was  about  three  cycles  into  its  life  where  the  system  had 
essentially  been  stabilized.  Io  re f  is  obtained  to  be  155  mA/g 
using  non-linear  parameter  estimation  [24,27]  and  the  value 
of  Z)j  is  found  to  be  1.25  x  10-9  cm2/s.  Note  that,  when  an 
initial  guess  of  10  or  1000  mA/g  is  assumed  for  /0,ren  it 
converges  to  155  mA/g.  When  an  initial  guess  of  10  1 1  cm2/ 
s  is  assumed,  D®  converges  to  1.25  x  10  9  cm2/s.  However, 
when  an  initial  guess  of  10~8  cm2/s  is  assumed  converges 
to  10  8  cm2/s.  Similarly,  when  an  initial  guess  10  7  cm2/s  is 
assumed,  D j  converges  to  1 0  7  cm2/s.  This  means  that  at  C/ 


8  discharge,  there  is  no  diffusion  limitation  and  the 
diffusion  coefficient  can  only  be  predicted  to  be  >  1.25  x 
10-9  cm2/s. 

The  simulated  discharge  curves  using  the  CDC,  VDC  and 
VDAC  model  at  the  rates  of  C/8  and  4C  are  displayed  in 
Fig.  4.  Note  that,  the  differences  in  the  predicted  discharge 
curves  by  the  VDC  model  and  by  the  CDC  model  are  only 
observed  at  the  discharge  rate  of  4C  and  above  and  that  even 
at  6C  rates,  the  approximate  solution  given  in  Eq.  (17) 
(VDCA)  matches  exactly  with  the  numerical  solution  of 
VDC  model.  Therefore,  for  galvanostatic  charge-discharge 
process  of  5  pm  particle,  under  any  rate  less  than  2C,  the 
simple  CDC  model  can  replace  the  VDC  model  and  for  rates 


Fig.  5.  The  effect  of  exchange  current  density  on  the  discharge  performance  of  graphite. 
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Fig.  6.  The  effect  of  diffusion  coefficient  on  the  discharge  performance  of  graphite. 


up  to  6C  approximate  solution  given  in  Eq.  (17)  can  be  used 
without  any  loss  of  accuracy.  For  particles  of  radius  10  pm, 
CDC  model  is  good  only  till  discharge  rates  of  C/2.  This  is 
true  because,  the  error  involved  in  the  CDC  model  is  based 
upon  '/y,  which  increases  four  times,  when  the  particle  size  is 
doubled.  And  hence,  the  CDC  model  is  good  only  for 
discharge  rates  less  than  2C/4  =  C/2.  Similarly,  for  a 
particle  size  of  10  pm,  VDCA  can  be  used  only  for  discharge 
rates  less  than  1.5C.  Note,  is  determined  to  be  0.0024- 
0.076  for  the  discharge  rate  ranging  from  C/8  to  4C  for  a 
5  pm  particle.  For  higher  rates  one  can  obtain  approximate 
and  accurate  solutions  using  higher  order  polynomial  profile 
approximations  [25,26].  Note  that,  at  higher  discharge  rates 
or  for  higher  particle  radius  CDC  model  cannot  be  used  and 
VDC  model  has  to  be  used  as  seen  for  a  4C  discharge  in 
Fig.  4b. 

Fig.  5  shows  the  discharge  potential  profiles  predicted  for 
the  Ni-composite  coated  graphite  at  different  reference 
exchange  current  densities.  We  observe  that,  the  lithium 
deintercalation  potential  decreases  with  increasing  exchange 
currents. 

The  simulated  lithium  discharge  potential  profiles  (for  a 
constant  C/8  discharge  rate  and  cut  off  potential)  with 
different  diffusion  coefficients  are  presented  in  Fig.  6  using 
VDC  model.  When  the  diffusion  coefficient  decreases  from 
1.25  x  10-9  to 2.25  x  10“u  cm2/s,  the  discharge  capacity 
decreases  from  316  to  306  mA  h/g.  The  decrease  of  the 
diffusion  coefficient  from  2.25  x  10“nto2.25x  10“ 12  cm2/ 
s  leads  to  a  large  decrease  of  the  capacity  from  306  to 
222  mA  h/g.  The  experimental  operating  condition  (Fig.  3) 
is  not  diffusion  limited.  One  can  only  conclude  that  the 
diffusion  coefficient  is  greater  than  1.25  x  10  9  cm2/s. 


5.  Conclusions 

A  simple  theoretical  model  is  presented  to  simulate  the 
galvanostatic  discharge  behavior  of  the  Ni-composite  gra¬ 
phite  electrode.  The  discharge  profiles  predicted  by  using  a 
constant  diffusion  coefficient  and  by  a  varied  diffusion 
coefficient  are  compared  in  this  paper.  The  results  show 
that,  the  VDC  model  can  be  simplified  to  the  CDC  model  for 
rates  less  than  2C  for  the  given  5  pm  particle.  Furthermore, 
for  rates  up  to  6C,  an  approximate  analytical  solution  is 
provided  for  the  VDC  model.  For  discharge  rates  less  than 
2C,  the  process  is  not  diffusion  limited  and  hence  diffusion 
coefficient  is  predicted  to  be  >1.25  x  ICC9  cm2/s  and  the 
exchange  current  is  predicted  to  be  155  mA/g. 
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